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Abstract 

We define natural probability measures on cycle-rooted spanning forests 
(CRSFs) on graphs embedded on a surface with a Riemannian metric. These 
measures arise from the Laplacian determinant and depend on the choice of a 
unitary connection on the tangent bundle to the surface. 

We show that, for a sequence of graphs (Q n ) conformally approximating the 
surface, the measures on CRSFs of Q n converge and give a limiting probabil- 
ity measure on finite multicurves (finite collections of pairwise disjoint simple 
closed curves) on the surface, independent of the approximating sequence. 

Wilson's algorithm for generating spanning trees on a graph generalizes to a 
cycle-popping algorithm for generating CRSFs for a general family of weights 
on the cycles. We use this to sample the above measures. The sampling 
algorithm, which relates these measures to the loop-erased random walk, is 
also used to prove tightness of the sequence of measures, a key step in the 
proof of their convergence. 

We set the framework for the study of these probability measures and their 
scaling limits and state some of their properties. 
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1 Introduction 



A cycle-rooted spanning forest (CRSF) on a graph Q is a subgraph each of whose 
components contains a unique cycle, or equivalently, contains as many vertices as 
edges, see Figure [TJ A cycle-rooted spanning tree (CRST) is a connected CRSF. 
Natural probability measures on CRSFs arising from the determinant of the graph 
Laplacian were introduced in [8]: the probability of a CRSF is proportional to the 
product of its edge weights times a product over its cycles of a certain function of 
the cycle, depending on the monodromy of a discrete C*- or SL2(C)-connection. 

The interest of these measures is that they give to a cycle a weight which is 
a function of its shape (homotopy type or enclosed curvature in the topological or 
geometrical settings, respectively, as described below). 




Figure 1: A CRSF on the 100 x 100 square grid; each connected component is in a 
different color and the cycles are in bold. 

We give a "cycle-popping" algorithm (Theorem [T]) for rapid exact sampling from 
these and other measures, generalizing the well-known cycle-popping algorithm of 
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Wilson [TT] for generating uniform spanning trees. One simply runs Wilson's algo- 
rithm, and when a cycle is created, flip a coin (with bias depending on the cycle 
weight) to decide whether to keep it or not. 

We use this sampling algorithm to sample approximations of scaling limits of 
CRSFs on surfaces in two different settings: a topological and a geometrical one, as 
follows. By surface we will mean here an oriented smooth surface with a Rieman- 
nian metric. By approximation of these, we will mean a sequence of graphs (Q n ) 
geodesically embedded on a surface E and conformally approximating E in a sense 
defined below (essentially, the simple random walk on Q n converges to Brownian 
motion on E). We endow the space of multiloops with a natural topology and show 
the weak convergence of probability measures on discrete multiloops to probability 
measures on continuous multiloops as follows. 

In the topological setting, we consider only the conformal class of the metric. 
An incompressible CRSF is a CRSF with no contractible cycles. Let (Q n ) be an 
approximating sequence and /j,f nc be the uniform measure on incompressible CRSFs 
of Q n . We show (Theorem 15) that the cycle process P™ nc of the /x" nc -random CRSF 
on Q n converges to a random loop process F inc on E, independent of the sequence Q n . 
This only depends on the conformal class of E, in the following sense. Let z±, . . . , Zk 
be distinct points of E, and S > be small; then for any isotopy class of sets of 
pairwise disjoint simple closed loops {71, . . . , 7 m } of E \ {.8,5(2:1) U • • • U B$(zk)}, the 
probability that a random incompressible CRSF on Q n has m cycles, and these are 
isotopic to the ji (which we call a cylindrical event), has a probability converging as 
n — > 00 to a limit independent of the approximating sequence Q n . This refines the 
result of [8j which shows that the distribution of the homotopy classes of the cycles 
in tti(E) has a conformally invariant limit. 

In the geometrical setting, we take into account the metric on E, and in particular 
its curvature. Let (Q n ) be a sequence of finite graphs conformally approximating E. 
Associated to this data is a discrete connection on a complex line bundle over Q n 
arising from the Levi-Civita connection on the tangent bundle on E. It is defined 
up to gauge equivalence by the property that the monodromy around a closed loop 
is e td where 9 is the enclosed curvature. From this connection $£c we construct 
a natural probability measure /i^c 011 CRSFs on Q n : each CRSF has a probability 
proportional to the product of its edge weights times the product over its cycles of 
2 — 2 cos 8, where 9 is the curvature enclosed. The corresponding loop process P^c is 
shown (Theorem 17) to converge to a probability measure Flc on multicurves on E. 

When the surface E is contractible, we define another measure which is in 
some sense more natural. This measure is a limit of the CRSF connection measures 
with curvature e ie9 when e — > 0, a limit which was introduced in [SI. This yields 
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a measure on CRSTs (CRSFs with one component) with weight proportional to 8 2 
where 9 is the enclosed curvature. We show (Theorem 17) that the loop measure 
P2c° conver g es to a measure Plc° on simple closed curves on S. 

In all three cases, our generalization of Wilson's algorithm is used to show tight- 
ness of the sequence of measures (Section |4.3 ). Convergence follows from the conver- 
gence of the measures on a determining class, the above-mentioned cylindrical events 
(Lemma 13) in the flat case and, in addition, convergence of polygonal approxima- 
tions in the curved case. 

We give samples from the measures ulc and (j,lc f° r the round sphere (Figures 
2 [3|, a saddle surface (Figure [5]) and a compact disk in the Poincare plane (Figures 
5), and for the measure /i inc on a flat torus (Figure [6J and planar domains (Figures [7] 
and [9]). For fx LC these are conditional samples, conditioned on having only loops 
with area (curvature) bounded by 7r/2; the above sampling algorithm does not work 
without this condition (see, however, [6] where it is shown how to sample from any 
determinantal process with Hermitian kernel, of which fiic is one). 

The paper is organized as follows. In Section [2] we introduce the sampling algo- 
rithm and prove its correctness. In Section [3] we introduce the probability measures 
on CRSFs on graphs on surfaces and show how they are exactly sampled by the 
algorithm. In Section [4] we show that the probability measures on loops that these 
induce converge to probability measures on the space of multiloops of the surface. 
Section [5] enumerates some of the properties of the measures on loops considered in 
the paper. The paper concludes with a list of open questions in Section |6l 



2 A general sampling algorithm 

An oriented CRSF is a CRSF in which each cycle has a chosen orientation. A 
measure on CRSFs induces a measure on oriented CRSFs by giving each cycle an 
independent 1/2 — 1/2-chosen orientation, and a measure on oriented CRSFs induces 
one on CRSFs by forgetting the orientation. 

Let Q = (V, E) be a finite graph with vertex and edge sets V and E, respec- 
tively, and c : E — >■ 1R >0 a positive function on the edges which we call the con- 
ductance. Let a be a function which assigns to each oriented simple loop 7 in Q 
a positive weight 0(7) G [0,1]. These functions c, a define a probability measure 
H = Hc,a on oriented CRSFs, giving an oriented CRSF Y a probability proportional 
to rieer c ( e ) llcycies 7 cr a {l)- We describe an algorithm to sample an oriented CRSF 
according to the measure /i. 

We note that this sampling algorithm requires a G [0, 1]; it will not work without 
modification for larger a. In the special case where a = 1 and c = 1, the algorithm 
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samples according to the uniform measure on oriented CRSFs. 

Let us describe a cycle-popping procedure, named P[w, T], which takes as argu- 
ments w a vertex and T an oriented subgraph of Q not containing w, and outputs 
another oriented subgraph of Q containing V and w. The procedure is the following: 
start at vertex w and perform a simple random walk (with each step proportional to 
the conductances) until it first reaches a vertex v which either belongs to T or is the 
first self-intersection of its path; 

• If v is in T, then replace T by the union of F and the oriented path just traced 
by the random walk. 

• If v is the first self-intersection, let 7 denote the oriented cycle thus obtained, 
and sample a {0, l}-Bernoulli random variable with success probability 0(7); 

— If the outcome is 1, then replace T by the union of T and the oriented 
path just traced by the random walk. 

— If the outcome is 0, erase the cycle that was just closed and continue to 
perform the random walk from v until it reaches T or self-intersects, in 
which case repeat the above instructions. 

The algorithm, called A, is then the following: start with T empty and w an 
arbitrary vertex in V, perform P[w, V] until T contains all vertices of Q, and output 

r. 

Theorem 1. If a (7) > for some 7 then the algorithm A terminates and its output 
is an oriented CRSF, sampled according to the measure \i. 

Remark 2. Note that if in the above definition of A one starts with V equal to S, a 
distinguished set of vertices in Q, then the algorithm samples an oriented "essential" 



CRSF with Dirichlet boundary conditions on S, see 18] or section 5.1 below for a 
definition. 

Proof. As in Wilson's proof, we construct an equivalent description, denoted A', of 
the algorithm A. 

Let us define A' in the following way. Consider, over each vertex v G V, an infinite 
sequence of i.i.d random variables = (x[ v \ x!f \ . . . ), each distributed as a 
random neighbor of v according to the conductance measure, that is, for each i > 1 
and each neighbor w of v, we have 

p (x^ = ^ - c{vw) 
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We represent as an infinite stack of cards, with only X{ v ' being visible at the 
top of the stack. 

We draw an edge from each vertex v to the neighbor shown on the top of the 
stack X^; the oriented graph thus seen is an oriented CRSF (with possible loops of 
length 2). This is our initial CRSF. We now describe a step by step random popping 
algorithm of the cycles. Note that at each step, the graph that we see remains an 
oriented CRSF. Here is the algorithm: For each cycle 7 encountered in the current 
CRSF, pop it with probability proportional to 1 — a (7); when a cycle is popped 
off, the top card on the stacks for each of its vertices is discarded. When a cycle is 
"kept" , its cards are fixed and can no longer be removed. The algorithm stops once 
the cycles that remain have been all "kept" in a Bernoulli trial. It is easy to see that 
the order in which the cycles are popped is not relevant. 

Note that this algorithm terminates (as long as there is at least one cycle with 
positive probability) because eventually, with probability 1, all the cycles present at 
one step will have been previously kept. 

Since the card of the stacks are distributed as the steps of a conductance-biased 
random walk, we see that algorithm A' has the same output in distribution as al- 
gorithm A. In order to compute the output distribution of A we will therefore use 
algorithm A'. 

Let us compute the probability that a given oriented CRSF T is obtained as 
an output of algorithm A'. Let 71,..., 7& be the cycles of T. The CRSF T is 
obtained as an output if and only if there exists a finite sequence of oriented cycles 
Ci, • • • , C m such that these cycles are popped, and after removing them, the cards 
that appear correspond to T, and there are k successful trials for Bernoullis with 
success probability «(7i). 

By independence of the cards in the stacks, the last CRSF considered is inde- 
pendent of the cycles that were popped. Therefore, for any oriented CRSF T, we 
have 

P(r) = P(r I pop C)P(pop C) 

C = {Cl,...,C m } 

= J^P(pop C, and T occurs underneath and is kept) 
c 

k 

= ^p(po P c)np(e)n a ^) 

c eer i=i 

= (^p(popc))np(e)n«( 7J ) 

V c J eer i=i 
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which we see is proportional to the weight of T. 



□ 



To sample a non-oriented CRSF according to a measure which assigns a CRSF 
T a weight proportional to rieer c ( e ) l~l 7 cr a (7)> where the product is over non- 
oriented cycles 7, and a is a function invariant under orientation, it suffices to have 
a G [0,2], perform algorithm A for the measure fi c ,a/2, and forget the orientation 
in the resulting oriented CRSF. In particular, we obtain the uniform measure on 
non-oriented CRSFs with the choice c = a = 1 . 

There is a variant of the previous algorithm to sample an oriented CRSF according 
to measure /i CjCf conditional on having a single loop: multiply all the loop weights by 
a small constant e. Then perform A; if e is small there will typically be a single loop 
(if not, start over). 

Let N be the total number of vertices of the graph. The running time of the 
algorithm is bounded by the time to obtain the first loop (which is bounded by 0(N 2 ) 
if a > 0(1/N 2 )) plus the running time of Wilson's algorithm, that is 0(iV(log N) 2 ) 
(Wilson's algorithm has a running time bounded by the cover time [TT] which is linear 
up to a logarithmic correction PQ). The running time is at least linear. Extreme cases 
correspond to extreme values of a: for a = 1 (uniform measure on CRSFs), the 
running time is linear; for a close to zero (like in the conditional measure described 
in the previous paragraph), the running time is large, at least 0(1/ sup a). 

3 Natural probability measures on CRSFs 

The most natural probability measure on CRSFs on a finite graph is the uniform mea- 
sure. There are however other natural probability measures that can be constructed 
from connections on line bundles and that are meaningful for graphs embedded in 
surfaces. 

3.1 Connections 

Let Q = (V, E) be a finite graph. A complex line bundle on Q is a copy C v of 
C associated to each vertex v G V. The total space of the bundle is the direct 
sum W = (& v( z V A unitary connection $ on W is the data consisting of, for each 
oriented edge e = vv', a unitary complex linear map (p vv i : C v — » Cv such that <p v r v = 
(p~ v , (equivalently, we assign to each oriented edge e a unit-modulus complex number 
ip e such that <£>_ e = (p e ). We say that two connections $' are gauge equivalent if 
there exist unit-modulus complex numbers ip v such that ipv'¥>w' = f'vv'^v, that is, 
is obtained from $ by changing the basis of each space C„ by a rotation. Let 
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c : E — y M >0 be a conductance function. We let A$ be the associated Laplacian 
acting on / e W denned, for each vertex v, by 

where the sum is over all neighbors v' of v. 

When Q is geodesically embedded in a surface £ (that is, embedded in such a way 
that edges are geodesic segments), there is a natural connection $ = $v on arising 
from a connection V on the tangent bundle TE: we define for each vertex v the line 
C v to be the tangent plane to S, with its natural complex structure coming from the 
metric. The V-parallel transport along edge e of Q defines the parallel transport <p e . 

3.2 Laplacian determinant and measures 

Theorem 3 ([51 IS])- For a graph with connection $ on a complex line bundle we 
have 

det(A*)= II c ( e ) II ( 2 - 2cos #)' 

CRSFs edges cycles 

where e %e is the monodromy of the connection around the cycle, for any choice of its 
orientation. 

Associated to $ is a probability measure /x$ on CRSFs, where the probability of 
a CRSF is proportional to Hedges c ( e ) ILyciest 2 - 2 cos 0). 

Another measure, introduced in [H], is obtained by passing to the limit t — ^ 
in a connection with weights ip e = e lt9e . This yields a measure on CRSTs (one- 
component CRSFs) which assigns a weight 8 2 to each CRST with cycle 7, where e ld 
is the monodromy around 7 (for any choice of its orientation) of the initial connection 
at t — 1. 



3.3 Flat connections 

Suppose that Q is geodesically embedded on a surface S with flat connection V on 
its tangent bundle. Let $ = $v be the associated connection on Q. Since V is flat, 
it has no curvature around contractible cycles. So the associated measure /i$ gives 
zero weight to contractible cycles and thus is supported on incompressible CRSFs. 

Let [ii nc be the uniform measure on incompressible CRSFs on Q. The measure 
//$ has density H 7 cr ( 2 — 2 cos #7) with respect to /i inc . 

Although /i inc cannot itself be written as a connection-measure //$ for some flat 



connection V, we can use the /i$ to study /i inc , see Lemma 13 below. 
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3.4 Graphs embedded on a curved surface 
3.4.1 The Levi-Civita measure 

Suppose that Q is geodesically embedded on a curved surface S. Take V to be the 
Levi-Civita connection associated to a metric g on S. Define fi^c to be the associated 
probability measure. It gives a CRSF a probability proportional to H 7 cr (2 — 2 cos # 7 ) 
where 9 is the curvature enclosed by 7. 



3.4.2 The CRST measure 

As discussed above, there is another measure Hlc° we can associate to this situation, 
but only when X is contractible. It is supported on the CRSTs of Q. Let $ = 
{e l8e } e( zE be the parallel transports on Q defined from V, and let $ t = {e Me } e( zE] 
these are well defined by contractibility of E. Let Ulc be the limit as t — > of the 
measures In Hlc°i eacn CRST T has a weight proportional to 6^, the square of 
the curvature enclosed by the unique cycle 7 of T. 

Let Z LC o = X^crsts 67 be the partition function of [lie - By Theorem |3j we have 

Z LC o = lim t~ 2 det A$ t . 

t— >o 

The partition function can also be computed as follows. Let Kg be the number of 



spanning trees of Q and T the transfer impedance function (see section 4.1.1 for the 
definition of transfer impedance). 

Theorem 4 ([7] Lemma 8). Let E + be an orientation of the edges. We have 
Zlc° = Kg X I #e - X G ^e'T(e, e) 

e&E+ \ e'£E+ 



3.4.3 CRST and LERW 

The measure \ilc° has a density with respect to a weighted "bridge" measure on 
loop-erased random walk excursions on the graph. Let Plc° be the measure on 
simple loops in Q obtained from fi^c by ignoring the branches of the CRST. Define 
Plerw t° be the law of a loop-erased random walk (LERW) from x to y where e = xy 
is a uniform random oriented edge in Q. For a simple loop 7, define Plerw(7) to 
be the probability that (1) the uniform random edge xy is an edge of 7, and (2) the 
LERW from x to y is the remainder of 7. 
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Lemma 5. For any cycle 7 in Q, we have 



Plco(7) = ^^rrPLERwM, 

^LC° |7l 



where \E\ is the number of edges of Q . 
Proof. We have 



V> LC o(j)Z LC o = e i 

jTD7 

I {spanning trees which contain j\ e}\ 



f 

7 

CRSTD7 

ItI 



eS7 

1 ^ r |^|K e P LE Rw(7). 



□ 



3.5 Exact sampling 

The measures /ij nc , /ilc°; an d /^lc can be sampled using our generalization of Wilson's 
algorithm as follows. 

The measure fi inc is sampled by using a function a which assigns a loop weight 
if it is contractible and 1 otherwise. 

For Hlc i we se t a (7) = for small e. For small enough e there will typically 
be only one loop (if there is more than one, start over). 

We can sample fi^c only in the case where the absolute value of the curvature 9 
enclosed by any curve doesn't exceed tt/2. Indeed, in that case, we will always have 
2 — 2 cos 6 G [0, 2] which is necessary to sample. 

See Figures |2j [3| |4j [5j [6] which are obtained by using fine conformal approximations 
to the underlying surfaces (only the cycles of the CRSFs are drawn). Figures [3] and 
[5] are samples conditional on enclosing curvature less than z/2. 

In order to sample the unconditional measures, one can use a general algorithm 
for determinantal processes with Hermitian kernels [B] which is slower. Figures [7] 
and [9] show samples of Hi nc on multiply-connected planar domains. 

4 Scaling limits for graphs on surfaces 

The measures fi inc , \ilc-, Hlc induce measures on sets of cycles on Q, by forgetting 
the rest of the CRSF. We use notation F inc , F LC , Plc° respectively for these cycle 
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Figure 2: A /i iC o-random CRSF on the sphere with its round metric 




Figure 3: A /z^c-random CRSF on the sphere with two components 
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Figure 4: A /^c- ran dom CRSF on on the saddle surface 
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Figure 5: A /x^c-random CRSF in a hyperbolic ball of radius 2 centered at (0, 1) 
the hyperbolic plane 
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Figure 6: A /ij„ c -random CRSF on the flat torus (which is obtained from the unit 
square by identifying opposite sides) 

measures. 

4.1 Conformal Approximation 

Let (Gn) be a sequence of (edge-weighted) graphs geodesically embedded in E with 
mesh size (longest edge length) going to zero. 

There are a number of equivalent definitions of the notion of conformal approxi- 
mation of £ by the sequence (G n )- Perhaps the easiest is to say that (conductance- 
weighted) random walk on G n converges to the Brownian motion on E, up to repa- 
rameterization. Another, more computationally useful, approach uses the transfer 
impedance. 
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Figure 7: A /i inc -random CRSF on a punctured disk conditioned on having a partic- 
ular homotopy type 
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4.1.1 Transfer impedance 

For any n let T n (e, e') be the transfer impedance of two oriented edges e = xy and 
e' = x'y' in Q n . This is defined as the (algebraic) current through e! when one unit 
of current enters at e_ and leaves at e + (the endpoints of e). In terms of the Green 
function G n one has 

T n (e, e') = G n (e+, e' + ) - G n (e + , e'_) - G n (e_, e' + ) + G n (e_, e'_). 

The function T™ is a discrete one-form on Q n (function on oriented edges which 
changes sign under change of orientation). We say that conformal approximation 
holds if when e±, e 2 are not within o(l) of each other, 

T n (e 1: e 2 ) = (A.A-g D (z lt z 2 )^ i(e 1 )l{e 9 ) + o{l{e x )t{e 2 )) , (1) 

where go is the continuous Green function, ^(e») is the edge length, and the notation 
tttQd represents the directional derivative in the direction of the edge applied to 
z { . 

If we sum for e% varying along a path from vertex a to 6, and e 2 on 7 Cj ^ 
from c to cf, we find that 

E G'"( e r,e 2 )-G"(er,e+)-G"(e+,e 2 ) + G"(e+,e+) = 

</r>(a, c) - £d(&, c) - 5- D (a, d) + g- (6, d) + o(l). (2) 

This quantity represents the change in potential from c to d when one unit of current 
enters at a and exits at b. As a function F(c) of c, this is the unique harmonic 
function with appropriate logarithmic singularities at a and b, and which is zero 
at d. In particular the level curves of F are equipotentials and, taking equipotentials 
near a and b (which are close to circles), one can thus compute the resistance between 
a small circle around a and a small circle around b. In this way, using convergence 
of the transfer impedance, one can show that the main notions of potential theory 
including the Poisson kernel, Cauchy kernel, holomorphic functions, etc. all behave 
well under conformal approximation. 

As an example of a conformally approximating sequence, we can pull back a fine 
square grid in C under a conformal map from a domain in C to E. More generally the 
(almost-)equally-spaced real and imaginary curves of a quadratic differential <p(z)dz 2 
on E naturally give a graph structure on E which is (almost-) conformally equivalent 
to a square grid graph except for conical singularities (with angles multiples of 2n) 
at the zeros of <j). 
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4.2 The probability space of simple multiloops 

For a graph Q n , let P n be one of the measures on loops discussed above. We see 
the probability measures P™ as living on the same probability space that we now 
describe. 

A simple multiloop in X is a finite union of simple loops (ignoring parameteriza- 
tion) with disjoint images, that is, an injective map from the union of k copies of the 
unit circle T to S, for some k > 1 (and modulo reparameterizations). The space of 
single simple loops is a metric space: the distance is defined by 

d(f,g) = inf sup \\f(t)-g(a(t))\\ 

01 tes 1 

where the infimum is over all reparameterizations a. In other words two loops are 
close if they can be parameterized so that their images are close for all parameter 
values t. One can extend this distance to the case of multiloops by taking the infimum 
over all permutations of loops, of the pairwise distances (and defining the distance 
between and fl^i for k ^ k' to be infinite). With this distance Q is a topological 
metric space. It is a disjoint union Q = U^L^k where Q k consists of multiloops with 
k components. 

This space is not complete: it is easy to construct Cauchy sequences that shrink 
to a point or to non self-avoiding loops. However it is separable: take all finite 
multiloops lying on fine lattice approximations of the surface. This is a countable 
family of loops which is dense in fi. 

The set of cycles of a CRSF on Q n defines an element of f2. 

A finite lamination on a surface is an isotopy class of a simple multiloop. For 
any points x±, . . . ,x m G X, and small 5 > 0, let Bi be a ball around Xi for some 
radius less than 5, and consider the finite laminations in E \ {Bi U ■ • ■ U B m }. Any 
simple multiloop 7 which avoids the balls Bi defines a lamination [7]. For any of 
these laminations L we consider the event 

E Bl ,..., Bm ;L = {7 e n I [7] = L} . 

We call these sets cylindrical events and consider the a-field B on Q generated 
by these events. 

Lemma 6. B contains the Borel sets in Q. 

Proof. We just prove this for one loop, that is, for fii. The proof is easily extended 
to the general case. 

Let c be a smooth simple closed curve in E and for some small e > let U e (c) be 
its e-neighborhood. Consider the event E that the random curve 7 maps into U e (c), 
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winding once around the annulus with, say, the positive orientation. These types of 
events generate the Borel sets. 

Let 21,22, • • • be a sequence of points dense in the boundary of U e (c). Let Bi 
be a ball around X{ of radius 5/i. Let Eg- n be the cylinder event that the curve 
7 in E \ {£>! U • • • U B m ) separates the points {21, . . . , x n } on the two boundary 
components of U e (c), that is, is consistent with 7 winding once around U E (c). The 
event E is contained in the intersection over n of the Eg. n . In fact we have E = 
r\ n E 0]n = Us^o H n Es-n- any continuous simple loop which separates the points lies in 
the interior of the annulus and winds once around. This can be seen as follows. First 
of all, any curve in n n E . n is contained in U e (c): otherwise, we could find some Xj 
lying on the wrong side of the curve since the family (x n ) n >i is dense in the boundary 
which would yield a contradiction. Second, the curve cannot be contractible since 
this would contradict the fact that it separates the points. Hence it winds once 
around the hole of the annulus. Its orientation is necessarily positive. □ 



4.3 Asymptotic size of loops and tightness of the measures 
4.3.1 Macroscopic loops 

For fi nc , the loops necessarily are macroscopic since they are noncontractible. In the 
curved case for the measures Flc and Flc , we show that there are, with positive 
probability, macroscopic loops. 

Theorem 7. With positive probability Flc o>nd Flc contain a macroscopic loop, 
that is, for sufficiently small e > the probability that there is a loop with diameter 
> e does not tend to zero with n. 

Proof. We show this for F L c°- Absolute continuity between it and Flc implies that 
it will hold for Flc as well. 

Consider first the n x n grid H n scaled by 1/n to the unit square [0,1] 2 . Let 
z\ 7^ z 2 G (0, l) 2 and let fx, f 2 be faces of the grid close to z\, z 2 respectively. In [7J, 
it is shown that the probability F(f 1 ,f 2 ) that fi,f 2 are enclosed in the cycle of a 
uniform CRST of H n is const /n 2 . Let E be the event that fx, f 2 are enclosed in the 
cycle. On this event, the area of the enclosing cycle is with high probability > en 2 
for some e > (see below, where it is shown that the loops are locally SLE2). Then 

for constants C, C. Here the denominator of the central inequality follows from the 
result of [7] that the second moment of the area of a uniform CRST is of order n 2 
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times a constant. The right-hand side of (J3J) is bounded by below by a positive 
constant independently of n. 

A similar argument holds for any graph Q n conformally approximating a surface 
E, and near a point where the curvature of the metric is nonzero, since the arguments 
of [7] which we used (in both the numerator and denominator above) rely only on 
the transfer impedance. □ 

We are not able to show that all loops are macroscopic with probability 1 (al- 
though we believe it is true). This appears to rely on a more delicate computation. 



4.3.2 Microscopic loops 

The following theorem shows that loops of P£c do not shrink to points as n — > oo. 
Theorem 8. Any subsequential limit ofF^ as n — > oo is supported on Q, that is 
limsupP^c (there is a loop of area < a) = o(a) . (4) 

n 

Proof. Let us argue by contradiction. If we suppose that Q is false, it means that 
there exists p > and arbitrarily small values of 5 > such that for n large enough 
(depending on 5), we have 

¥ n LC (there is a loop with area > p . 

Hence, using the Markov property and the cycle-popping algorithm, we show that 
there exists < q < p and arbitrarily small values of 5 > such that for n large 
enough (depending on 5), we have 

¥ n LC (there are 1/5 loops with area < 5 2 ) > q 1/s . 

However, by Lemma [27] below, there exists < £ < q such that for any k and n large 
enough, we have 

P2 C (there are > k loops) < £ fe . 

By taking k = 1/5 we obtain that for arbitrarily small values of 5 > 0, there is a 
large enough n such that 

< q 1/s < ¥ n LC (there are 1/5 loops with area < 5 2 ) < £ 1/s . 

This yields a contradiction since the right-hand side of the last equation tends to 
zero faster than the left-hand side when 5 — > 0. □ 
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4.3.3 Resampling and tightness 

We will show tightness of the sequence of measures P™ nc , P^o and P^c*- This will 
yield the existence of subsequential limits. 

We show in fact that the scaling limits of the macroscopic loops are absolutely 
continuous with respect to SLE 2 . Let 7 be a loop in a P2c°~>^£.c~ or ^tnc' 1 'andom 
CRST, and a, b 6 7 distinct points on it. Let j[a, b] be the part of 7 counterclockwise 
between a and b, and 7 [6, a] the complementary part. If we erase j[a, b], we can define 
a new loop 7' by taking a LERW from a to b in the domain defined by £ U 7 [6, a], 
with wired boundary conditions on 7 [6, a], and conditioning on the LERW to end 
at b. The union of this LERW and 7 [6, a] is the simple closed curve 7'. By the 
sampling algorithm, this curve 7' is absolutely continuous with respect to 7 [a, b], 
with probability ratio given by the ratio of the cycle weights. If a and b are close 
to each other, the LERW from a to b will with high probability not exit a small 
ball around [a, b]. Hence the cycle weight of 7' will be close to that of 7. Wilson's 
algorithm thus shows that this is a fair sample of P conditioned on 7 [b, a}. Thus the 
LERW from a to b with the appropriate boundary conditions is absolutely continuous 
with respect to 7 [a, b], with a bound independent of mesh size 1/n. 

This proves that in the scaling limit, the loops are locally absolutely continuous 
with respect to the scaling limit of LERW, which was shown in [10J to be SLE2. 

In particular the scaling limit is supported on simple curves (recall that SLE 2 
curves are simple, and note that this is a local property). 

We also need to check that for P" nc and P^c the number of loops doesn't blow up 



as n — > 00. This is a consequence of Lemma [27] below. 

Theorem 9. The sequence (P^ c ) n >g is tight on f2. 

The above argument also shows that the sequence (P£ c o)n>o is tight, provided 
we allow for the possibility that the curve shrinks to a point. We need to add to 
Qi a copy of £ whose points represent the constant maps of S 1 to that point. The 
metric naturally extends to this augmented space so that the limit of a sequence 
of curves shrinking to a point is the constant curve at that point. Let be this 
augmented space. 

Theorem 10. The sequence (P2c°) n >o ^ s ^S^t on f£J\ 

It may be that (P2c°) n >o * s tight on fii but we have not proved that its curves 
do not shrink to points as n — > 00. 

By Theorem [8] above, (P£c) n >o * s tight on f2, that is, its curves do not shrink to 
points. 

Theorem 11. The sequence (P2c)n>o ^ s ^idht on f2. 
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4.4 Probabilities of cylindrical events 

4.4.1 Flat connections 

Let Xi,...,Xk € £. Let M. be the space of flat SU 2 -connections modulo gauge 
transformations on E\ {xi, ■ ■ ■ ,Xk}- Such a flat connection is determined uniquely 
by a homomorphism from 7Ti of the surface into SLV 

The fundamental group 7Ti(E\{xi, . . . , Xk}) is a free group F m on m = 2g+k—l+b 
letters, where g is the genus and b the number of boundary components of E. Thus a 
flat connection is determined by m arbitrary elements of SU2, one for each generator 
of 7Ti. Let v be the canonical measure on Ai: the image of Haar measure on SU™. 

4.4.2 Trivalent graph 

A useful device, see [J], is to define a trivalent graph H in E (unrelated to Q) with 
a single puncture of X in each face, so that if is a deformation retract of S. We can 
thus think of £ as a ribbon graph structure on H . 

Recall that a finite lamination L of £ is an isotopy class of a finite number of 
simple curves. A lamination L retracts to a multicurve (with "multiplicity" ) on H; 
L is determined by, for each edge of H, a nonnegative integer giving the number of 
strands of a minimal representative of L retracted onto that edge. These integers 
satisfy the conditions that at each vertex of H the sum of the three integers is even 
and the three integers satisfy the triangle inequality: see Figure [8j Moreover any set 
of integers satisfying these two conditions arises from a unique lamination. 



Figure 8: Strands of the lamination at a vertex of the trivalent graph; the integers 
i, j, k satisfy i + j + k = mod 2 and \i — j\ < k < i + j . 




k 
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We define a partial order on laminations: L < V if the integers on edges of H 
associated to L are all less than or equal to those of V . We define the complexity 
n(L) of L to be the sum of these integers. 

An SU2-connection on H determines a flat connection on S. After gauge transfor- 
mation one can take the SU2-connection to be the identity on all edges of a spanning 
tree of H. 

4.4.3 Integrals over the space of flat connections 

Given a finite weighted graph Q embedded in £, and $GM define 

= e n c ( e ) n ( 2 - Tr K)) , (5) 

CRSFs edges cycles 7 

where cj 7 is the monodromy of the connection <3> around the cycle 7, for either choice 
of its orientation (and any starting vertex). The function Z is a real- valued function 
on M. the space of flat connections modulo gauge. 
We denote 

z ° = e n c ( e )> 

inc. CRSFs edges 

the partition function for all incompressible CRSFs, without cycle weight. (This is 
not the same as Z(Id), which is zero.) 

For a flat connection $ we may rewrite ^ as 

Z($) = ]Tx L T L ($), 

L 

where the sum is over finite laminations L, where Xl is the conductance- weighted 
sum of CRSFs whose cycles are isotopic to L, and 

T L ($) = ;Q(2-Tr(u; 7 )) 

76 L 

is a real- valued function on A4. 

Fock and Goncharov [1] proved that, seen as real- valued functions on Ai, the 
functions Sl = Yl^eL Tr(w 7 ) are linearly independent and generate the vector space 
of polynomial functions on A4, when L runs through all finite laminations. Hence 
the functions T L = Yl^eL (2 ~ Tr(w 7 )) are also linearly independent, generate the 
same vector space and, when the bases are ordered by increasing number of cycles, 
the change of basis matrix M$ y T is an invertible infinite triangular matrix. 
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Choose an ordering of the Tl consistent with the partial order on laminations 
defined above. Let P = {Pl($)} be the Gram- Schmidt orthonormalization of the 
T = {Ti($)} with respect to this ordering and with respect to the inner product 
(/)#) = Jm fgdv. Let A = (Al,l') be the infinite lower-triangular matrix such that 
P = AT. 

Recall the linear operator A$ on the total space C 2 ' y ' of the C 2 -bundle on Q. 
Theorem 12 ([8]). We have 

Z($) 2 = det(A#). 

One can extract the coefficients of any desired lamination L as follows. 
Lemma 13. For any cylindrical event El, we have 



L'>L 

This sum is finite for any finite graph. 
Proof. The probability of E L is X L /Z . Write 

Z($) = J^XlTl = X*T($) = X'A _1 P($) . 

L 

Since P is orthonormal, we have J M ptpdv = Id. Hence, 

/ ZY l Kdv = I X*A- 1 P($)P*($)Arfz/ 
Jm Jm 

X'A" 1 (^J P($)P*($) du^j A = X* . 



Hence, 



Each entry is the integral 



X= / ZA f Pdu. 

I M 



X L = J2 A L',L [ Z($)P L ,($)du. 
r;->r. J M 



L>>L 

Dividing by Z we obtain the result. □ 
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4.5 Convergence in the flat case 

We first consider $ to be a flat connection. Associated to this is a measure /xj 
on incompressible CRSFs. Since /i| has a density (independent of the graph) with 
respect to /i" nc , it suffices to show that this latter converges. 

The main tool is the following convergence result. Let x±, . . . ,Xk be points of E 
and Bj a small ball around Xj. Let $ be a flat connection on E \ {Bi U ■ • • U -B&}. 

Theorem 14 ([9]). There exists a function F e L 2 (A4) depending only on the 
conformal type of the surface E \ {Bi U ■ ■ ■ U B^} such that 



The proof follows from Theorem [12] and the convergence of discrete harmonic 
functions to continuous ones. 

Theorem 15. Let Q n be a sequence of graphs with mesh size going to zero and 
conformally approximating a non-simply connected surface E. The measures P n = 
P™ nc which are uniform on incompressible CRSFs of Q n converge as n — )■ oo ; and the 
limit is a conformally invariant measure F inc supported on sets of pairwise disjoint, 
noncontractible simple loops on £. 

In [5] the homotopy classes on E of the noncontractible loops were shown to have 
a conformally invariant limit distribution. 

Proof. Take points z±, . . . , in E, take 5 > small, and for each i let Bi be the ball 
of radius 5 around Zj. Let Q"% = Q n \ {Bi U ■ ■ ■ U B^} and the associated measure 
on multiloops of incompressible CRSFs on Q n whose loops stay in Qg. A loop on 
£ \ {B\ U • ■ ■ U Bk} is said to be peripheral if it is isotopic to one of the boundary 
curves dB{. Laminations on Q*% may contain peripheral loops, but we are interested 
in those without peripheral loops. 

Up to errors uniform in n and tending to zero with 5, 

F n {E L ) = F n B (E L | no peripheral loops) . 

This follows from the sampling algorithm since removing one or more very small 
holes does not change the distribution of the LERW. 

It remains to show that lim^oo P^(^) exists and depends only on the conformal 
type of the domain E \ {B\ U ■■• U the conditioning multiplies by a non-zero 
constant which is also a conformal invariant quantity and independent of L. 
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Consider L a finite lamination in E \ {B\ U • • • U B^}, which has no peripheral 



curves. By Lemma 13 the probability P^(El) is given by a sum over V > L of 



integrals over Ai of Z(Q)/Zq times a function P/y(<3>) independent of n. 

By Theorem 14 the integrand Z(Q)/Zq converges. Moreover, it is bounded inde- 
pendently of n as follows. We write 

Z*/Z = J2 x l/Z T l . (6) 

L 

First T L < e ° (n(L)) because 2 - Tr( w) is uniformly bounded by a constant over Ai. 



Secondly, by Lemma 27 below, X% < e~ cn ^ Zq for arbitrarily small c > 0. Since 
there are at most n(K) M laminations L with complexity n(L) = n(K), where M is 
the number of edges of H, the sum ^ is bounded by a converging series. Hence by 
bounded convergence, for each V the integral J M ^^-P^dv converges. 
We need to show that the sum over V converges. We write 

(£X K T K ,P V ) = (J2 x kTk,Pv) (7) 

K K>U 

since (T Kl P v ) = unless K > V . We also have \(T K ,P L ,)\ < e oi -< K ^ by the 
Cauchy-Schwarz inequality. Using the above bounds the sum Q is bounded by 
e - c ' n ( L )Z for some arbitrarily small d > and summing over L' now gives a conver- 
gent sum. 

We now use a classical convergence argument [2] . We have shown that the prob- 
abilities of any cylindrical event converge. The cylindrical events form a family of 
sets which is stable under finite intersection. Since it generates the a-field B it is a 
determining class, that is if two probability measures on (O, B) coincide on all the 
cylindrical events, then they are equal. 

Since we furthermore have tightness by Theorem |9j the sequence of probability 
measures P n admits subsequential limits by Prokhorov's theorem. More precisely, 
for any subsequence, there exists a subsubsequence which converges. Since its value 
on the cylindrical events is known, there is only one possible limit. Let us call it F inc . 
Now, since Q is a metric space, this implies that the sequence /j^ nc converges weakly 
to F inc . □ 

As a corollary, we obtain the convergence of the measures /xj for any flat connec- 
tion $. 

Corollary 16. There exists a probability measure P$ on (£1,13) such that 
in the sense of weak convergence. 
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Proof. It suffices to show the convergence of this sequence of measures on finite 
intersections of balls in Q of small radius because this is a determining class for B. 
For any curve 7, there is a small radius r > such that its tubular r-neighborhood 
retracts onto 7. Any curve in this r-neighborhood and winding once around is isotopic 
to 7. On any such neighborhood the density YI^cl (2 — w 7 — u~ : ) is constant, hence 



the convergence follows by Theorem 15 □ 



4.6 Curved case 

The measures ^lc^lc converge in the following sense. 

Theorem 17. There exist probability measures Flc'Plc on (0, jB) and (£l*,B) re- 
spectively such that for any sequence (Q n ) of graphs, geodesically embedded on E ; with 
mesh size going to zero as n — )■ 00 and conformally approximating E ; the sequences 
of probability measures ¥^ c and P^c conver 9 e weakly towards respectively F LC and 

Proof. Let us approximate E by a polygonal surface E e , that is, with a surface which 
is flat except for conical singularities. A standard way to do this is to take a fine 
triangulation of the surface (with triangles of diameter at most e and whose angles 
are bounded from below), and replace each triangle with the Euclidean triangle with 
the same edge lengths. As e — > the conformal structure on E e converges to that of 
E. (Indeed, there is a homeomorphism from E e to E which is (1 + o(l))-biLipschitz.) 

We furthermore suppose that the conical singularities are separated at distance 
6(e) ->• with s = 0(6). 

Any graph embedded on E or E e can be embedded on E e or E with small dis- 
tortion; furthermore a graph conformally approximating E e will have image on E 
conformally approximating E, and vice versa. 

Let Zx, . . . , Zk be the vertices of E e . The Levi-Civita connection on E e is a flat 
connection on E e \ {z\, . . . , Zk} and approximates the Levi-Civita connection on E, 
in the sense that the curvature enclosed by any loop is close for both connections 
(for small e). Restricting to Q n , this shows that P2c e is close to P^c (since cylinder 
events have close measure). 



By Corollary 16, if we fix E £ and take Q n embedded on E £ and conformally 
approximating it as n — > 00, the measures P2c e converge as n — > 00 to a limit Plc.e- 
Similarly for P^c-o £ - 

For any e > the probabilities of the cylindrical events (away from the singu- 



larities) are determined by the function F of Lemma 14 These are expressed as 



integrals of the Green function on paths avoiding the singularities. The addition of a 
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singularity modifies the Green's function, and hence the probability, by a negligible 
amount. Since the singularities are at distance 5 » e it follows that the probability 
measures fic,e and Pz,c°,e eac h form a Cauchy sequence. 

By a diagonal argument, we can take n — > oo and e — Y and conclude that P^c 
and P^c converge weakly as n — > oo. □ 

From experimental simulations, the probability of getting two or more loops in 
Ulc on the round sphere is on the order of one percent. Hence measures ulc and 
fi^c are close for the total variation distance. 

The measure F LC o is a limit of Flc when the metric is scaled by a factor t — > 0. 
Hence it is not obvious that Theorem [8] implies that Plc° also is supported on 
macroscopic loops. However, we conjecture it to be true. 

Conjecture 18. F LC o supported on fl\ that is 

Plc° (the area of the loop is zero) = . 



4.7 A comment and two applications 

The convergence result of probability measures /i Cja on CRSFs is actually more gen- 
eral and can be adapted for a wide range of functions a on the cycles, not necessarily 
coming from connections. This is due to the fact that the crucial convergence argu- 



ment is made for the uniform measure on incompressible CRSFs (Theorem 15). We 
now state two applications of the previous theorems. 



As an application of Theorem [T7J we obtain a positive answer to Conjecture 1 
of [?]. Let Q n be a rectilinear approximation of a planar Jordan domain D (see [7] 
for a definition). Let A denote the combinatorial area (number of faces) of the cycle 
of a cycle-rooted spanning tree of Q n . Let 6 = A/n 2 be the Euclidean area of the 
cycle. We denote by E" nif the expectation with respect to the uniform measure on 
CRSTs on Q n . 

Corollary 19. For k > 2 there exists ctfc(-D) > such that 

K mf (A k )=a k (D)n^(l + o(l)). 
Proof. Let k > 2. We have 



= n 2fe "Xco (6 k - 2 ) E n uni{ (A 2 ) 

= n 2k - 2 C(D)\D\E LC o (9 k - 2 ) (1 + o(l)) 
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where the last equality follows from Theorem 8 of [7j (here C(D) is a constant) and 
the weak convergence of P^c ^° P^c - Since 6 is bounded and for P^c is with 
positive probability non-zero by Theorem |7[ the limit E £C o (# fc ~ 2 ) is a positive real. 
The corollary is proved by taking a k (D) = C(D)\D\E LC o (6 k ' 2 ). □ 



As an application of Theorem 15, let us state the following corollary 



Corollary 20. Consider a uniform spanning forest on an annulus-graph wired on 
its boundary. Then the simple closed curve separating the two connected components 
has a conformally invariant limit which is given by the measure Pj nc conditional on 
having only one component. 

Proof. This follows from the fact that the dual of a wired essential forest on the 
annulus is a uniform incompressible CRST on the annulus with free boundary con- 
ditions. The measure is thus given by /i" nc conditional on having one loop. The 



convergence follows from Theorem 15 □ 



Figure [9] shows a sample of this interface between the two tree components of a 
wired uniform spanning forest on the annulus. 



5 Properties of the measures 

In this section, we mention a few properties of the measures on CRSFs on surfaces 
we have been considering. This could help give a characterization of these measures. 

5.1 Markov property 

CRSFs on surfaces satisfy the following spatial Markov property. Consider a graph 
G embedded in a compact oriented surface E of positive genus. Let a : — > M >0 
be any positive weight function on the cycles of Q. 

Let {7} = {71, . . . , 7fc} be a family of cycles in Q which separate £ in a number 
of connected components Si, ... , E r , with r < k. For i — 1, . . . , r, denote by Gi the 
intersection of Q and the closure of Sj (i.e. Sj along with its boundary) and by dGi 
the boundary cycles. 

An essential CRSF with Dirichlet boundary conditions on a graph with boundary 
is a spanning subgraph which is the union of disjoint unicycles not touching the 
boundary and a union of trees rooted on the boundary. 

Let r be the random CRSF on Q associated to the connection $. For each i, 
denote by Tj the random essential CRSF on the line bundle Gi, with boundary dGi, 
associated to $j obtained by restriction of $ on Gi- 
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Figure 9: A uniform incompressible CRST on an annulus in the square grid 



Proposition 21. Conditional on {7} C T, the random CRSF F is equal in distri- 
bution to 

{7} u[ =1 r, . 

Proof. The key argument is that, because of the assumption on {7}, a CRSF r on Q 
which contains {7} is the union {7} U r i=1 Tj of {7} with essential CRSFs Tj on each 
one of the connected component S, of S \ {7}. The proposition then follows directly 
from the cycle popping algorithm. □ 



5.2 Restriction property 

Let Di C D be two planar Jordan domains. Let Q be a finite graph approximation 
of D. Let $ be a connection on a line bundle over Q. We denote det A D (<3>) the line 
bundle Laplacian on Q with connection $. For any subset S of the set of vertices 
of Q, we denote by det A^*($) the line bundle Laplacian with Dirichlet boundary 
conditions on S, see [8]. 
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Lemma 22. For any finite set of simple non-intersecting curves {7} C D\, we have 

/ff ({7}) _ ( jet A^($) \ / det A g ($) \ 
^ ({7}) ^detAg } ($) J 1 ^detAf 7} (d>) ) " 

5.3 Stochastic domination 

Recall from Section [2] the notation /i CjQ , for measures that assign a CRSF V a proba- 
bility proportional to rieer c ( e ) l~l 7 cr a (7)- 

Lemma 23. Lei Si C 5*2 be two subgraphs of Q . Let Pi and P2 6e essential CRSF 
measures fi C)Ct with Dirichlet boundary conditions on Si and S2, respectively. For any 
curve 7 G Q \ S2, we have 

Pi(7)>P 2 (7)- 

Proof. This follows from the definition of algorithm A described to sample these 
measures since it is less likely for LERW, started from a vertex outside 5*2, to create 
loop 7 before it hits 5* 2 than before it hits S\. □ 

Note that in the case the measures come from a line bundle connection $, this 
implies 

detAg lU7 ($) > det Ag 2U7 ($) 
detA Sl ($) ~ detA Sa ($) ' 

which is a non trivial potential theoretic consideration (which can be translated in 
terms of Dirichlet-to-Neumann map). 

5.4 Sub- exponential tail for the distribution of the number 
of loops 

Let P n be one of the two sequences of measures P™ nc or P^c*- We show that the number 
of loops has tail decreasing faster than exponentially. This involves a computation 
from two specific examples, one for Flc and one for Pm C . 

5.4.1 Wired boundary for Flc 

Let Q be a cylindrical square grid graph with height n and circumference n, obtained 
from a (n + 1) x n square grid by identifying the left and right boundaries. We wire 
the upper and lower boundaries, that is, add an extra vertex connected to all points 
along the upper and lower boundary. 
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Let $ be a line bundle connection on Q with monodromy q = e 1 ^ around each 
square face, where = c/n 2 for a constant c. Thus the total curvature is c. 

We can realize $ by putting parallel transport q y from (x, y) to (x + l,y) (indices 
cyclic) and the identity on vertical edges. We compute the partition function Z(Q) 
of the measure on essential CRSFs of Q n (CRSFs with wired boundary). 

Theorem 24. = 1 + 0(c 2 ). 

Corollary 25. With probability 1 — 0(c 2 ) there are no loops. 

Proof of Theorem. The matrix of the Dirichlet Laplacian A$ commutes with hori- 
zontal translation (rotation of the cylinder). The action of A$ on the e^-eigenspace 
(here 9 = ^ for k = 0, . . . , n — 1) for horizontal translation is 



A*(0) 



/4-2cos# -1 \ 

-1 4-2cos(fl + </>) -1 



-1 ••. -1 

\ -1 4-2cos(9 + (n-l)(f))J 

For small we have A$ = A/ + S for a small diagonal perturbation S, so 

det A$ = det A 7 det(J + SAj 1 ). 

Keeping terms of order up to 2 gives 

Z§ det A$ 
^7 ~~ det A/ ~~ 

= II f 1 + E ^) a fJ#) + E a,:.^)a,;.^! - A,:^)A,>); 

Here the product is over 9 = 2-rrk/n, k — 0, . . . , n — 1, and 

S^.(0) = 2j(f) sin + f(f) 2 cos + O(0 3 ). 

The entries of Aj x (9) are written in terms of Chebyshev polynomials. It is 
convenient to define a > by 4 — 2 cos 6> = 2 cosh a. Then for i > j we have 

! = sinh(jq) sinh((n + 1 - i)a) 
I,hj sinho;sinh(n + l)a 

(Since the matrix is symmetric this defines the case i < j as well.) 
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If a is not small, AjJ -(0) ~ e ^ i ' a /2 sinha. For a small we can use instead the 
bonnd Aj^e) < 

Expanding out the product over 9 we find that the first-order term in <fr vanishes 
because sin 9 is odd. We are left with 



where 



- 7 ^ = l + <j ) 2 (A 1 + A 2 + A 3 ) + 0( ( j ) 3 



Ai = y^^. 2 sinh(jq) sinh((n + 1 - j)a) ^ Q 
sinha sinhfn + I) a 

9 j=l V ; 

A * = EE 4 ^ sin2 ^ [ A 7h( e ) A Jid d ) - a tU 6 ) a iXM > 

e j<k 



and 



^3= E E 4 ^ sin ^ sin ^ A ?]j A w 

0i<9 2 j,k 

e j,k 

Since = c/n 2 , it remains to show that A\ + A 2 + A3 = 0(n 4 ). 
Summing A 2 and A 3 gives 

a 2 + a 3 = - e E 4 ^ sin2 £, ( A 7W a - ^ E E 4 ^' 2 sin2 ^ A ili) 2 - ( 9 ) 

9 j<k e j 

The second sum in (|9| is 0(n 4 ). The first is 0(n 4 ) if we restrict to \9\ > e for a 
constant e. Similarly the sum ([8]) is 0(n 4 ) under this restriction. It now remains to 
consider the case of \9\ < e. 

Consider first A\. The sum over j can be computed explicitly as a power series; 
the leading term is of order n 3 and is 

n 3 cos9 cosh(n + l)a 
6 sinha sinh(ra + l)a 

When 9 is not 0(l/n), we have ^^("+1)" ~ ^> anc ^ s i n ' 1 « = 9 + 0{9) 2 so the leading 
term is 
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When 9 = 2nj /n and j is of order 1, cos M n + 1 ) Q: j s Q f constant order and so the terms 
for small 9 contribute a negligible amount to the final sum. 

Similarly the sum over j in ^ can be computed explicitly (we resorted to using 
a computer algebra program); the leading term is of order n 3 and is 

n 3 sin 2 9 _ 3 / 1 
~Y (e 2a - l)sinha ~ " \Q9 

for 9 not on the order of 1/n. 

The leading terms in A 1 and A 2 + A 3 thus cancel and the sum is n 3 ^2 e O(l) = 
0{n A ). 

Thus 

^ = l + 4<V) = l + 0(c 2 ). 
Zji n 

□ 

A similar computation holds for a cylinder with aspect ratio m/n — > t. By 
universality, the results holds for any graph conformally approximating the curved 
annulus. If the curvature is not constant the result still holds since the measures for 
the constant curvature and variable curvature case are absolutely continuous with 
bounded Radon-Nikodym derivative between them. 

A rectangle with wired boundary can be obtained from an annulus with wired 
boundary by wiring the vertices along a zipper connecting the two boundary com- 
ponents. Since wiring increases the probability of having no loops, we can conclude 
that for a rectangle with total curvature c, the probability of having no loops is also 
1 — 0(c 2 ). By universality and conformal invariance this holds for any disk and 
similarly for any planar Riemann surface. 

5.4.2 Wired boundary for P inc 

As above let Q be a cylindrical square grid graph with height m and circumference n, 
obtained from a (n + 1) x m square grid by identifying the left and right boundaries. 

Let $ be a flat line bundle on Q with monodromy z around the boundary. We 
can realize $ by putting parallel transport u from (x, y) to (x + 1, y) (indices cyclic) 
and the identity on vertical edges, where u n = z. 

We compute the partition function Z(z) of the measure on incompressible 
essential CRSFs of Q (CRSFs with wired boundary). 
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Theorem 26. As n, m — > oo with n/m — > r, the 'probability generating junction for 
the number of loops tends to 



Pr(X) = n 



q 3 + q 



2 + X 



3=1 



qj + q-j - 1 



where q = e 7TT . 

In particular the probability that there are no loops tends to a nonzero constant. 
The analogous computation with free boundary was done in [8]. 

Proof. By [8 J we have 

Z(z) = Y[ Ch n (X + 2) - z - 1/z 
x 

where Ch n is a variant of the Chebyshev polynomial defined by Ch n (a + 1/a) = 
a n + a~ n and A runs over the eigenvalues of the line graph Laplacian with wired 
endpoints, that is, eigenvalues of the m x m matrix 

\ 



( 2 
\ 



-1 

1 '•• -1 
-1 2 



These eigenvalues are 2 — 2 cos for j — 1 



, re- 



setting X = 2 — z — 1/z in Z(z), the power of X records the number of loops in 
the configuration. We have 

Z(z) = J]Ch n (A + 2)-2 + X, 

A 

and the probability generating function for the number of loops is 

Ch n (4-2cos^r) -2 + X 



n 



m+l ' 



Ch n (4-2cos 



TTJ_ 



(10) 



j=l " """"m+l' 

When n is large, Ch n (w) is exponentially large in n for real w unless w is near 2. 
Thus the X dependence of (10) arises from j near 0. Letting n, m — > oo with 
n/m — > t this expression is asymptotic to (with q = e 1TT ) 

^ q j + q- j -2 + X 



n 

3=1 



q j + q~ 



1 



□ 
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5.4.3 Loop number lemma 



Lemma 27. For any < £ < 1, there exists N and K such that for all n > N and 
for all k > K, we have 

P n (there are at least k loops) < £ fe . 

Proof. Each loop created during the performance of the cycle-popping algorithm 
either disconnects the surface or decreases the rank of the first homology. By the 
Markov property, the law of the conditional CRSF is obtained by independently 
sampling in each of the connected components. 

For P2c> suppose that we have created k closed curves. Each of the complemen- 
tary components is either planar or nonplanar. Each new loop found has a positive 
probability of being macroscopic, that is, either reduces the rank of the first homol- 
ogy or removes definite area from both resulting components, by Theorem [7j Thus 
for k sufficiently large, the curvature enclosed c on either side will be small, hence 
with definite probability 1 — p — 1 — 0(c 2 ) the two sides will contain no further loops 



(by Corollary 25 above). By taking k large enough, and mesh size 1/n small enough, 
the curvature enclosed c can be taken as small as needed such that p < £. 

Thus the probability of an extra loop eventually decays exponentially with rate 
less than £. 



A similar argument works for Pj nc using Theorem instead of Corollary 25 . □ 



In fact we conjecture that the decay of the probability of the number of loops 
decays faster than exponentially since each new loop subdivides a region into regions 
each having curvature smaller by a definite fraction. This would imply a decay of 
order £ fc jk\. We conjecture an even stronger decay as follows. 

Conjecture 28. Let P be one of the measures F^c or F inc . There exists < £ < 1 
such that for all k > 1, we have 

P (there are at least k loops) < £ k . 



6 Questions 

1. Can the measures F inc , F LC o, and P^c* be defined directly instead of via limits of 
CRSF measures? For example via a stochastic differential equation, like a vari- 
ant of SLE 2 defined on Riemann surfaces subject to some potential depending 
on the metric? 
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On the round sphere for the measure (jllc -, can one use t ne connection Lapla- 
cian to say more about the shape of the cycle, as is done in [7j in the flat case? 



By Corollary 19 above, the expected area of a P iC o-random loop is jgrpjrpj > 



where a^(D)n is the leading order of the third moment of the combinato- 
rial area of the loop of a uniform cycle-rooted spanning tree on Q n . Can one 
compute a^{D)l Is there a number theoretic interpretation? 

3. What can be said about the Gaussian Free Field associated to the line bundle 
Laplacian? Are our loop models related to this GFF? In particular, for a choice 
of an infinite curvature, we expect to obtain loops at all scales. 

4. What is the right framework for the study of Plc? What is the distribution of 
the number of loops? 

5. Is there a probabilistic interpretation for the coefficients of the triangular ma- 
trix A defined in Section 14.4.31 ? 
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